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ABSTRACT 

It has long been thought that there is a connection between ultraluminous 
infrared galaxies (ULIRGs), quasars, and major mergers. Indeed, simulations 
show that major mergers are capable of triggering massive starbursts and quasars. 
However, observations by the Herschel Space Observatory suggest that, at least 
at high redshift, there may not always be a simple causal connection between 
ULIRGs and mergers. Here, we combine an evolving merger-triggered AGN 
luminosity function with a merger-triggered starburst model to calculate the 
maximum contribution of major mergers to the ULIRG population. We find 
that major mergers can account for the entire local population of ULIRGs hosting 
AGN and ~25% of the total local ULIRG luminosity density. By z ~ 1, major 
mergers can no longer account for the luminosity density of ULIRGs hosting AGN 
and contribute <12% of the total ULIRG luminosity density. This drop is likely 
due to high redshift galaxies being more gas rich and therefore able to achieve high 
star formation rates through secular evolution. Additionally, we find that major 
mergers can account for the local population of warm ULIRGs. This suggests 
that selecting high redshift warm ULIRGs will allow for the identification of high 
redshift merger-triggered ULIRGs. As major mergers are likely to trigger very 
highly obscured AGN, a significant fraction of the high redshift warm ULIRG 
population may host Compton thick AGN. 

Subject headings: galaxies: active — quasars: general — galaxies: Seyfert - 
galaxies: starburst — infrared: galaxies 



1. Introduction 



In the 1980s astronomers discovered a new class of infrared selected galaxies known as 
ultraluminous infrared galaxies (ULIRGs) and characterized by Ljpt > 10 12 L , where Ljpt 
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is the 8-1000 fim luminosity (e.g., iHouck et al.lll984j ; ISoifer et al.lll984j ). Another important 



topic during this time period was the study of the evolution of quasars (e.g.. lSchmidt &: Green 
19831 ) . a class of active galactic nuclei (AGN) where accretion onto the supermassive black 
hole at the center of a massive galaxy gives rise to Lx > 10 44 erg s _1 , where Lx is the 
2-10 keV luminosity. ULIRGs and quasars have similar bolometric luminosities (10 45 -10 46 
erg s -1 ) and optical observations suggest many ULIRGs have nuclear sources of non-thermal 



lomzmt 



g radiation and disturbed morphologies ( panders et al.lll988al ). Thus, ISanders et al. 



( 1l988al ) suggested that when two gas rich galaxies merge, gas and dust will fall into the 
nucleus of the resulting galaxy, triggering a massive starburst and a quasar. 

Thirty years later, the connection between ULIRGs, AGN, and major mergers is still 
an area of active research. Recent simulations show that, indeed, gas rich major merg- 



ers are capable of triggering large starbu r sts and bright AGN (e .g., iHopkins et al.l l2006al ; 
Younger et aDl2009l : IHopkins et allhoiObl : iNaravanan et al.lboioh . By looking for ULIRGs 
with strong X-ray emission or a power-law spectra in the Spitzer Space Telescope IRAC 



bands, studies have shown that AGN are common in ULIRGs fe.g.. IHopkins et al. 



2006a 



Armus et a"Dl2007l : lYounger et~al1l2009l ; bonlev et aDboioi IHopkins et alihoiOri INaravanan et al 
2010l ). Furthermore, both the fraction of ULIRGs that host AGN and the fraction of ULIRGs 
whose bolometric lu minosities are dominated by AGN emission, appear to increase strongly 



with l uminosity (e.g.. IVeilleux et al 



20061 : iDonlev et al.l l2010l : iNardini et al. 



2002 



Pearson! 120051 ; iGruppioni et al.ll2005l ; iBrand et al. 



20101 ). Morphological studies have shown that a 



significant fraction of ULIRGs have disturbed morphologies, sugge s ting the galaxy has re- 
cently undergone a merger or int eraction (e.g. Donley et al.l 1201 Or iKartaltepe et al.l 12010 



Nardini et al.ll2010l ). Additionally. IVeilleux et al. 



(120021 ) showed that the fraction of ULIRGs 



triggered by major mergers increases with luminosity. Thus, it is expected that a significant 
fraction of ULIRGs host AGN and were triggered by gas rich major mergers. However, these 
studies tend to focus on ULIRGs with z < 1. 

Far-infrared observations by the Herschel Space Telescope have opened a new window 
on the z > 1 ULIRG population. Interestingly, Herschel o bservations show t hat at high 
redshift major mergers are not necessary to trigger ULIRGs dSturm et al.ll2010l ). Analyzing 
Herschel observations of the Bootes field, Melbourne et al.l (120121 ) point out that <30% 
of optically-faint z ~ 2 ULIRGs show obvious signs of a recent merger. Deep Herschel 
observations of the Great Observatories Origins Deep Survey (GOODS) and the Cosmological 
Evolution Survey (COS MOS) fields find that most z > 1 ULIRGs are not in a starburst 



mode of star formation (jElbaz et al.ll201ll ; iRodighiero et al.ll201ll ). Instead, the increased 
gas fr action in high redshift galaxies allows normal secular star formation to power ULIRG s 



Daddi et al 



20081 : lElbaz et al.l 1201 ll : IRodighiero et al.l l201lt Melbourne et all 12012 ) . 



(e.g., 

Recent cosmological simulations confirm these observational results, finding that more than 
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half of hi gh redshift ULIRG s can be accounted for through mechanisms other than major 
mergers (INiemi et al.l 120121 ) . These results indicate that at high redshift the AGN-ULIRG 
connection may be quite different than the connection observed locally. 



Draper fc Ballantynd ( 120121 ) computed an AGN population model that constrains the 



space density and Eddington ratio evolution of AGN triggered by major mergers by consid- 
ering the hard X-ray luminosity function (HXLF), X-ray AGN number counts, the X-ray 
background, and the local mass density of supermassive black holes. Thus, a model of the 
major merger population can be computed by combi ning this descr i ption o f the evolving 
luminosity function of merger-triggered AGN with the iHopkins et al.l (l2010bl ) model for the 
time evolution of merger-triggered starbursts. Similarly, a model merger spectral energy 
distribution (SED) can be calculate d by combining AG N infrared spectra computed with 
the photoionization code Cloudy (IFerland et al.l Il998l ) with the iRieke et al.l (120091 ) star 
formation templates. This method is used here to determine the maximum contribution of 
mergers to the ULIRG population at z < 1.5. A ACDM cosmology is assumed with Ho = 70 
km s" 1 Mpc" 1 and fi A = 1.0 - Q M = 0.7. 



2. Calculations 



2.1. AGN Model 



d2012|) 



The evolving HXLF of major merger-triggered AGN presented by lDraper fc Ballantyne 
is used to determine the spa ce density and luminosity distribution of merger-triggered 



AGN. 



Draper fc Ballantynd (120121) combin ed an observationally motivated m ajor merger rate 



( Hopkins et al.lboiOa ; Treister et al. I boioh and a theoretical AGN light curve (IHopkins &: Hernquistl 



2009[) to determine the contribution of major merger-triggered AGN to the HXLF. IHopkins fc Hernquist 
( 120091 ) parametrize the AGN light curve as 



\(t) 



\t-t Q \ 



1/2 



-2//3 



(1) 



where X(t) is the AGN Eddington ratio at time t years after the AGN was triggered, £q = 
t§rf I (2/3 In 10), and t , i], and (3 are para meters describing the AGN life time, peak Eddington 
ratio, and light curve slope, respectively. iDraper fc Ballantynd ( 120121 ) find that the best fit 
to the observed HXLF, X-ray AGN number counts, X-ray background, and mass density of 
supermassive black holes is achieved when to — 2.5 x 10 8 yrs, rj = 2.5, and = 0.7 for AGN 
triggered by mergers. It is assumed that all mergers trigger AGN and therefore this model 
provides an upper limit to the contribution of major mergers to the AGN population. 
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Here, this merger-triggered AGN HXLF is converted t o an infrared lumino sity function 
using AGN SE Ds computed using Cl o udy v ersion C08.00 (IFerland et al.lll998l ) as described 
in Section 2 of iDraper &: Ballantynd (1201 ll ) . These SEDs cover the sub- mm to very hard 
X-ray wavelength regimes and incorporate the direct emission from the AGN, the diffuse 
emission radiated along the line of sight by the obscuring material aro und the AGN, and the 
emissi on reflected off the inner surface of the obscuring cloud. As in IDraper fc Ballantyne 
(120111 ) , the neutral hydrogen density of the clouds is assigned such that Compton thin clouds 
have nu = 10 4 cm -3 and Compton thick (CT) clouds have Uh = 10 6 cm -3 , in accordance 
with the observed densities of typical molecular clouds. The inner radius of the obscuring 
material is assumed to be ~10 pc. As we compare against observations of an ensemble of 
sources, instead of fitting individual sources, a simple AGN torus model is appropriate for 
this study. Moreover, as discussed in Section [31 the infrared emission of ULIRGs is dominated 
by star formation processes, thus the results presented here are not dependent on the torus 
model used. 



2.2. Starburst Model 



The evolving s tar formation rate, M* (t), of the merger-triggered starburst is determined 
using the model of iHopkins et al.l (l2010bl ). Thus, 



dt 
dM* 



£* ln(10) exp 



-M* 



M sh /U 



(2) 



where t* = 0.1 Gyr is the timescale of the merger-triggered starburst and M s i, is the total 
mass of stars born during the merger -triggered starburst. Once M* is determined, Ljr due 
to star formation is calculated using ( THopkins etldlhoiOrJ fl 



Lir = 1.1 x 10 



10 



IMoyr- 1 



(3) 



In order to combine the AGN and starburst models, the time delay between the trig- 
gering of the starburst and the triggering of the AGN must be considered. Observational 
studies and simulations show that A tends to peak ~100 Mvr after M* peaks (jDayies et al. 



20071 : ISchawinski et al.l 120071 . 120091 : iRodrfguez Zaurin et all 12010 : IWild et al.ll2010l : IHopkins 



1 Using the iKennicuttl (|1998f ) conversion factor, 5.8 xlO 9 L,0/(M0 yr 1 ), yields an infrared luminosity 
density ~10% lower. 
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20121 ). Therefore, the starburst is triggered such that the time delay between the peak star 
formation rate, M*! eak , and the peak AGN Eddington ratio, A peafc , At = 100 Myr. Scenarios 
where At = and 1 Gyr are also considered. Once the merger occurs, the AGN light curve 
of Equation [T] and the star formation rate evolution of Equation [2] determine the evolution 
of the system. 

Thus, the space density of major mergers at redshift z is set by the merger rate. The 
starburst model is triggered at z + Az, where Az is the appropriate change in redshift for At 
at z. The corresponding AGN is triggered at z. The evolution of the starburst is governed 
by Equation [2] and the AGN luminosity evolves according to Equation [IJ The infrared 
luminosity of the merg er remnant galaxy is then determined using the A GN SEDs calculated 
with Cloudy and the iRieke et al.l (120091 ) starburst SED templates. The iRieke et al.1 (120091 ) 
templates are based on the SEDs of local, pure star forming galaxies and include poly-cyclic 
aromatic hydrocarbon (PAH) emission features. Figure [1] shows an example SED at four 
different times, t, after the major merger: before M* reaches Mf eafc (upper left), M* « Mf afc 
(upper right), A ~ \P eak (lower left), and after A reaches \P eak (lower right). With the merger- 
triggered AGN HXLF, g?$ x / 'd(\og L x) , and merger SEDs set, the merger infrared luminosity 
function, d$i R /d(logL IR ), is calculated as 



d$ IR (L IR , z) = d$ x (L x ,z) d{\ogL x ) 
d(\ogL IR ) d(hgL x ) d(\ogL IR ) 



(4) 



where Lx and Lj R are computed using the combined AGN-starburst SEDs. The infrared 
luminosity density, is then computed as 



^merger (^) 



J IR~ 



d§ IR (L IR , z) 
d(\ogL 



IR) 



d(\og L IR ), 



(5) 



where L"f R 



, 3 and Lfg x 



10 



13.5 



L m . To calculate the ULIRG ^ 



10 8 L e 

increased to 10 12 L Q . The ^ merger calculated here is then compare d with the in 
nosity density of AGN and their host galaxies, ^agn, measur ed by 
and the total infrared luminosity density, ^ total, measured by iLe Floc'h et al. 



merger j 



J IR 



IS 



rared lumi- 



(120051). 



20101. l201lh 



T he major merger rate is set by parametrizing observations and simulations (IHopkins et al 



2010al ; iTreister et al.l 120101) and the AGN light curv e parameters are set by fitting observa- 
tions of the AGN HXLF (IDraper fc Ballantyndl2012l ). The evol ution of the merger-tri ggered 
starburst model is parametrized by fitting results of simulations (IHopkins et al.ll2010bl ). This 
leaves only M sb and At as free parameters. We explore 9.0 < logM sfe /M < 11 in steps of 
0.25, assuming all mergers result in a similar value of M s b. The values At = 0, 100 Myr, 



and 1 Gyr are considered. The maximum ULIRG ^merger 
investigated. 



agn and ^ merger I ^ total are then 
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Results 



In order to prevent over-predicting the maximum ^merger, we consider t h e 24 / jm number 
count of X-ray selected AGN following Equation 1 of iDraper fc Ballantynd ( 1201 ll ). If M s b > 
10 10 - 25 M the bright end of the 24 lira number count is over-predicted by a factor >2. These 
same models also over-predict the local ULIRG \E* total- Thus, 10 10 - 25 M is an upper limit on 
the average M s j for the population of major mergers. Figure [2] shows the number count for 
mergers with A t = 100 Myr and M s b = 10 10 M . As the observed number count is from the 
GOODS fields (ITreister et al.ll2006l ). it is not surprising that the model slightly over-predicts 
the bright end of the observed count since GOODS is a narrow field survey and likely misses 
bright, rare sources which are better accounted for by wide field surveys. 

Figure [3] shows the maximum ULIRG ^ merger /^agn and ULIRG ^ merger I ^ total that 
does not significantly over-predicting the 24 [xm number count. For the model shown in 
Figure El At = 100 Myr and M sb = 10 10 M . The top left frame of Figure [3] shows that 
mergers can account for the local ULIRG ^agn, however, as shown in the lower left frame, 
^ merger I ^ total < 0.26 for local ULIRGs. Interestingly, if we remove the starburst from our 
SEDS, we find that emission from AGN alone can contribute <20% of the local ULIRG 
^ agn- Thus, we confirm that the Lir of ULIRGs hosting AGN tends to be dominated by 
star formation. 



If M s b = 10 10 M , major mergers can account for the local ULIRG ^ agn for At = 
and 1 Gyr. Thus, the ability for mergers to account for the local AGN ULIRG population 
is not strongly dependent on At. The local ULIRG ^ merger/^ total reduces to ~0.20 for 
both At = and 1 Gyr. If At = 100 Myr and M s b < 10 9 ' 75 M , mergers cannot account 
for the local population of ULIRGs hosting AGN. Observations and simulations suggest 
that M s b is proportional to the galaxy st ellar mass, M*, such t hat Ms b (M*) = fM * with 
the fraction / on the order of, but <0.1 (IHopkins et al.l l2010bl ; IZavala et al.ll2012l ). If we 
assume a distrib ution of M^fMj defined by M sfe (M„) = 0.1M*, with the distribution of M* 
described by the iPerez- Gonzalez et al.l (120081 ) stellar mass function, the AGN number count 
and the local ULIRG ^ total are over-predicted by at least a factor of 2. However, if / = 0.05, 
the model predictions are in agreement with the AGN number count and the local ULIRG 
^ merger I ^ total = 0.28, only 0.02 higher than the simple calculation with a single M s b- If 
/ = 0.05, (M ab (Mj) = 10 10 - 05 M . Therefore, for At = 0-1 Gyr and M sb ~ 10 10 M , the 



local ULIRG * 



merger 



AGN 



1.0 and the local ULIRG 



merger 



total 



0.20-0.28. 
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Discussion 



Major mergers can account for at most a quarter of the local ULIRG \f* total- This 
suggests that a large fraction of local ULIRGs are triggered by mechanisms other than 
the coalescence of two massive gas rich galaxies, such as minor mergers, interactions, and 
secular processes. By z ~ 1, the ULIRG ^ merger agn < 1 and by z ~ 1.25, the ULIRG 



m.erger/^ AGN 



merger I 

S, 0-5. At z ~ 1, the ULIRG ^ merger/^ total < 0.12. Indeed, simulations by 
Hopkins et al.l (l2010bl ) predict that at z = 1, the ULIRG \I/ mer (, er /^totai ~ 0.08, in good 
agreement with the findings of this study. As major mergers are much more common at 
z > 1 than z < 1, this suggests that secular processes are even more important for triggering 
ULIRGs at high redshift. Herschel has provide d observational evi dence that major mergers 
are not necessary at z > 1 to trigger ULIRG s ( jSturm et al.l 120101 ); a finding that has been 
confirmed by simulations (jNiemi et al. 2012 ). 



At z > 1 the majority of ULIRGs appear to be scaled up versions of local normal 
star forming g alaxies. These high re dshift ULIRGs tend to have cooler far-infrared dust 
temperatures ([Melbourne et al.l 120121 ) and stronger PAH emission (lElbaz et al.l 120111 ) than 
local ULIRGs. F urthermore, b y com paring the 8 /xm flux to the total infrared flux of high 
redshift ULIRGs, lElbaz et al.l ( 120111 ) find that most z >l ULIRGs lie on t he infrared main 



sequen ce a nd are in a normal st ar forming mode of evolution. Similarly, iBournaud et al. 
fcpl2h and lKocevski et al.l koid \ find that at z ~ 0.7 and z ~ 2, AGN t end to be hosted by 



galaxies on the main sequence of star formation (see lElbaz et al.l 120 111 ). Observations and 
simul ations both point out that high redshift galaxies are more ga s rich than local galaxies 



(e.g., iTreister et al.ll2010l ; |Pi Matteo et al.ll2012l ; iNiemi et al.ll2012h . It is likely that galaxies 
with large reservoirs of gas and dust are capable of fueling ULIRGs without being triggered 
by a major merger. 



Elbaz et al.l (120111 ) do find that some h igh redshift ULIRG s are in a starburst mode 



possibly triggered by a merger. Moreover, lElbaz et al.l ( 120111 ) find evidence for a popu- 
lation of very highly obscured AGN embedded in these compact dusty starbursts. Other 
observational studies also fin d CT AGN (N H > 10 24 cm" 2 ) candidat e s in high red shift 



dusty starburst galaxies (e.g., iDonley et al.ll2010l ; iMelbourne et al.ll2012l ). iFabianl ( 119991 ) ex- 



plains that the gas and dust funneled into the central regions of the merger remnant galaxy 
will fuel a burst of star formation, rapid black hole accretion, and will obscure the result- 
ing AGN. Thus, a significant fraction of CT AGN are e xpected to be recently trig gered, 
likely by a major merger, and accreting very rapidly (e.g., iDraper fc Ballantynell2010l ). The 
results of this study are consistent with the merger-triggered CT AGN scenario, but, as 
^ merger/^ total < 0.12 for ULIRGs at z ~ 1, the contribution of these merger-triggered CT 
AGN to the z > 1 ULIRG population must be fairly small. The fraction of AGN that are 
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CT is hard to observationally constrain due to the severe obscuration that defines CT AGN 
and AGN population models predict a wide range for the CT fraction (see iBallantyne et al. 
201ll ). Increasing the fraction of AGN that are CT in this model by a factor of 1.5 increases 



the z ~ 1 ULIRG * merger /* total by -0.01. 

Local ULIRG SEDs are often divided into two groups, warm and cool, where warm 



d2006l) 



ULIRG SEDs are characterized by /^m/ ' fm^m > 0.2 (e.g., lAlonso-Herrero et al.l 120 06). 
where is the 25 /xm flux and /eo^m is the 60 jttm flux. According to lAlonso-Herrero et al. 
warm ULIRGs are 15-30% of the Bright Galaxy Survey sources (ISanders et al. 
The merger model used here does produce warm ULIRG SEDs with f^m/ ft 



1988bl ). The merger model used here does produce warm ULIRG SEDs with /25 Mm / feofim ^ 
0.5 at all redshift s; thus mergers can account for the local population of warm ULIRGs. 
Elbaz et al.l (120111 ) found that galaxies on the main sequence of normal star formation, in- 
cluding galaxies hosting AGN, tend to have cooler dust temperatures than star-bursting 
galaxies. Because starbursts triggered by major mergers are expected to be more compact 
than secular star formation, merger-triggered sta rbursts will be cha r acterized by higher dus t 
temperatures than normal star formation (e.g., lElbaz et al.l l2011t iMelbourne et al.ll2012l ). 
Indeed, observations of local ULIRGs show that compac t ULIRGs are more likely to host 
an AGN than less compact UL IRGs (jNardini et al.ll2010l ). Analyzing simulations of major 
mergers, lYounger et al.l ( 2009 ) find that warm ULIRGs are likely to be galaxies evolving 



from the star formation dominated merger phase to the AGN-starburst post-merger phase. 
Thus, by combining Herschel and Spitzer or Wide-field Infrared Survey Explorer (WISE) ob- 
servations to select warm ULIRGs, the population of ULIRGs hosting merger-triggered AGN 
can be identified. As discussed above, the population of ULIRGs hosting AGN triggered by 
major mergers is likely to include a significant fraction of CT AGN. Thus seeking out high 
redshift ULIRGs with warm SEDs will also lead to the identification of high redshift CT 
AGN. 

By combi ning the evolving lum i nosity function of AGN triggered by major mergers 
calculated by iDraper fc Ballantynd ( 120121 ) with the merger-triggered starburst model of 
Hopkins et al.l (j2010bl ). we computed an upper limit for the major merger contribution 
to the ULIRG population. Locally, major mergers can account for the observed pop- 
ulation of ULIRGs hosting AGN and ULIRGs with warm SEDs, but the local ULIRG 
^ merger/^ total ^ 0.26. By z ~ 1, major merger-triggered ULIRGs hosting AGN can no longer 
account for the population of ULIRGs observed to have AGN signatures. Indeed, the ULIRG 
* merger /* AGN < 0.50 by z ~ 1.25. Furthermore, at z ~ 1, the ULIRG * merger , < total < 0.12. 
Combining observations by Herschel and Spitzer or WISE to identify high redshift ULIRGs 
with warm SEDs is a good tool for identifying the population of ULIRGs hosting merger- 
triggered AGN, a large fraction of which are expected to be CT. 
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Fig. 1. — Infrared SEDs of an AGN-starburst galaxy with black hole mass M, = 10 8 M and 
obscuring column density Nh ~ 10 23 cm -2 , at various times after being triggered by a major 
merger. The dashed blues lines sh ow the AGN infrar ed SEDs computed with Cloudy and 
the red dot-dashed lines show the iRieke et al.l ( 120091 ) starburst SEDs. The solid black lines 
show the SED of the AGN-starburst galaxy. In the upper left plot, the AGN and starburst 
both contribute approximately half of Lir. The upper right plot shows the SED at M* w 
M^ eak , and thus the starburst dominates Lir. The lower left plot shows the SED at A ~ 
\P eak . The AGN dominates Ljr in both of the bottom plots. 
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Fig. 2. — 24 /im number counts for mergers with At = 100 Myr, M s b 
Lx > 10 42 erg s _1 . The blue lines show the ULIRG number count and the green lines show 
the LIRG (LiR > 10 11 L ) number count. The black lines show the 24 /im number count of 
all mergers. The solid lines show the number counts for major mergers and the dashed lines 
show the number counts for the merger-triggered AGN only. Data poi nts show the number 
count of X-ray selected AGN in the GOODS field fTTreister et aJbood ). 
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Fig. 3. — Fraction of ^agn and ^ total contributed by major mergers with At = 100 Myr and 
M s b = 10 10 M Q . The top row shows the fraction ^agn that can be accounted for by mergers 
for ULIRGs, LIRGs, and so urces with L j p > 10 8 L^ . The observed ^agn is determined by 
fitting the data points from iGoto et al.l (120101 120111 ). whic h are shown as t he cy an points. 
The bottom row shows the fraction of ^ total, as reported by lLe Floc'h et al. J2005L that can 
be accounted for by major mergers for ULIRGs, LIRGs, and sources with Ljr > 10 8 L . 



